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INTRODUCTION 


In  an  investigation  of  active  magnetic  detection,  a  basic  issue  is  the  decoupling  of  the 
source  field  from  the  detector.  One  concept  for  accomplishing  this  is  a  search  coil  gradiometer. 
TTie  source  coil  is  wrapped  around  the  center  of  a  permeable  core,  and  two  detection  coils  are 
wrapped  around  the  core,  one  near  each  end,  and  connected  in  series  opposition.  The  detection 
coils  and/or  the  source  coil  can  be  made  movable  along  the  axis  of  the  core,  in  order  to  achieve 
optimum  nulling  of  the  source-detector  coupling.  In  addition  to  this,  the  device  can  be  used  in  a 
discriminator  configuration,  in  which  the  source  fr^uency  is  rejected  in  the  detector  by  means  of 
a  feedback  circuit  with  some  time  constant  appropriate  to  the  application.  Of  course,  this  mode 
of  operation  requires  relative  motion  between  the  detector  and  the  scattering  object,  but  in  many 
applications  this  is  acceptable. 

At  practical  distances  from  a  scattering  object,  the  scattered  field  can  be  approximated  by  a 
linear  combination  of  a  uniform  field  and  a  uniform  gradient.  The  uniform  gradient  is  the 
quantity  to  be  measured  by  the  search  coil  gradiometer.  It  is  well  known  that  the  presence  of  a 
permeable  core  will  cause  magnetic  flux  concentration,  and  will  enhance  the  measurement  of  the 
magnetic  field.  The  question  to  be  addressed  here  is;  what  does  the  permeable  core  do  to  the 
field  gradient? 


THE  MODEL 


To  investigate  the  effects  of  a  permeable  core,  the  core  will  be  approximated  by  a 
permeable  prolate  spheroid.  The  spheroid  is  oriented  with  its  symmetry  axis  along  the  z-axis, 
and  a  uniform  field  plus  a  uniform  gradient  field  is  applied.  The  magnetic  potential  for  the 
external  field  is  then  given  by 


<I)‘  =  -H\x 


■  z^)  -  Gj'jxy  -  Gijxz  -  Gjsy  z  +  C 


(1) 


where  C  is  an  arbitrary  constant.  The  form  given  here  produces  a  unform  gradient  tensor  which 
is  symmetric  and  traceless,  consistent  with  the  requirements  V  •  H  =  0  and  V  x  H  =  0. 

To  solve  the  boundary  value  problem,  it  is  necessary  to  change  to  prolate  spheroidal 
coordinates  (^,T1,<|))  where  ^(1  <  ^  <  <»)  labels  a  family  of  confocal  prolate  spheroids, 

Ti(-1  <  T]  <  1)  labels  an  orthogonal  family  of  confocal  hyperboids,  and  (1)(0  <^<  2n)  labels 
angular  position  about  the  symmetry  axis.‘  With  the  foci  located  at  ±a  on  the  z-axis,  the 
connection  between  the  coordinate  systems  is  as  follows 


2a 


(2) 


1  "Static  and  Dynamic  Electricity",  by  William  R.  Smythe,  Third  Ediuon,  McGraw-Hill,  1968 


css  TM  567-91 


<t»=  fan  \x,y). 


where 


ri  =  yx^  +  y^  +  (2  -haf 


r^^yx^  +  y^  +  iz-af. 


The  inverse  relationships  are 


x=a'S(%^- 1)(1  -il^)cos(j) 


y=a‘\(^^-  1)(I  -Ti^)sin(|) 


z=a^Tl. 

The  general  solution  to  Laplace’s  equation  in  prolate  spheroidal  coordinates  is' 


0=  I  I[A,„cosm(j)-hB,„sinm(!)][C,X"(^)+D,„er(^)][£,„Pr(Ti)-HF,„Qr(Ti)]  (10) 

/ sOm  =  0 

where  P”  and  QP  are  the  associated  Legendre  functions  of  the  first  and  second  kind, 
respectively.  Note  that  Qn^)  is  regular  for  ^  ->  <»  while  Pn^)  is  regular  for  ^  1.  The  first 

few  functions  have  the  explicit  forms 

n°(^)=l  (11) 
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P‘(4)  =  3^V(^"-1) 
/’2(^)  =  3(4'-1) 


e;©=V?^5>n(|^)-^' 

0"®  =  j(34"-l)ln[|^]-f 

c;®=VF^fi{|^)-|f^' 


G2^(^)  =  (^^-1)  -In 


3,  r^+0.5^-3^^ 


To  put  the  external  applied  potential  into  prolate  spheroidal  coordinates  it  is  necessary  to 
use  (Eq.  7-9),  using  (Eq.  1 1-16)  for  guidance.  This  gives 


j:=fl/’,(4)/’,'(Tl)cos<j) 

(23) 

y  =aE,‘(4)E,‘(T|)sin(l) 

(24) 

z=aP”(^)P®(q) 

(25) 

xz=^/’2(4)F2'(q)cos(l> 

(26) 
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yz=jP2i^)PUr\)sin^  (27) 

and 

2 

xy  =^P^,(^)P^(r\)sin2^  (28) 

for  the  obvious  cases.  To  handle  the  terms  involving  squares  of  the  coordinate  values  it  is  more 
convienient  to  use  linear  combinations.  First 

ix^-z^)-(y^-z^)=x^-y^  =  ^Pl{%)Pl{y\)  cos  20.  (29) 

The  other  combination  is 

(x"-z')-(-Cy^-z")=x"-i-/-2z"  (30) 

=  fl\^^-l)(l-V)-2flW 

=  -y(34^-l)(3Ti^-l)-~ 

The  additive  constant  in  (Eq.  30)  can  be  absorbed  into  the  constant  in  (Eq.  1),  so  it  can  be 
dropped  without  loss  of  generality.  With  the  results  given  in  (Eq.  29-30),  the  following 
expressions  result 

2  2 

^'-z'  =  ^/'2(^)/’2(t1)cos20-^P°(^)F°(ti)  (31) 


and 

2  2 

y'  -  z'  =  -^Plik)Plmos  20  -  ^P”(^)F°(ti).  (32) 


The  external  potential  can  now  be  expressed  in  prolate  spheroidal  coordinates  as 
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l  =  2m  =  l 

s  X  [/4^'„cos//i<J)-(-5,'„sinm(j)]/""(^)P,'"(Ti) 

/  =  1  m  =  0 

(33) 

where 

^1,0  ~  ^2.0  ~  0 

(34) 

Al,  =  ~aH;  ,  ^,%=|^(C,',  +  Ca 

(35) 

A;,=-aH;  , 

(36) 

(37) 

K2=^(0^-g;,)  .  B«=^G,'a. 

(38) 

The  total  external  potential  has  the  form 

<D  =  0*-i-0‘’ 

(39) 

where  the  spheroid-induced  pan  0°  has  the  form 

0°=  i  I  [A,‘;„cosm<!)-»-S"„sin/n(l)]j2"(^)P"(Ti). 

/  =  1  m  =  0 

(40) 

The  interior  potential  0'  is  given  by 


l  =  2m  =  l 

<I>‘=  I  I  [/i;,„cosm(j)  +  B;^sinm(j)]P"(^)P"(il). 

i  =  1  m  sO 


(41) 


The  coefficients  are  determined  by  applying  the  boundary  conditions  of  continuity  of  the 
tangential  H  field  and  normal  B  field  at  the  surface  of  the  spheroid  where  ^  =  The  H  field  is 
given  by  H  =  -VO,  and  in  prolate  spheroidal  coordinates,  this  gives 
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and 


where  the  scale  functions  are 


and 


/i2dn 


U  1 


h,=  I  - - 

The  boundary  conditions  are 


and 


H^  +  H(  =  xHi\ 


(42) 

(42) 


(43) 


(44) 

(45) 


(46) 


(47) 


(48) 


where  x  = 


Upon  application  of  (Eq.  47)  after  deleting  common  factors  and  using  orthogonality  of  the 
trigonometric  functions,  there  results 
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AlMm  -Al„PT^)  =  -Al„PT(^)  (49) 

and 

The  remaining  equations  are  obtained  from  (Eq.  48),  using  orthogonality  of  the  trigonometric 
functions: 

and 

-  xBi^pni) = -Buni).  isi) 


Define  the  determinant  of  coefficients  D,  „  and  the  Wronskian  W, respectively,  as 


D,,„=pr&67(0-tpr&Qr& 


and 


Then  the  solutions  for  the  coefftcients  are  given  by 


Al„  = 


Bl.= 


and 


(53) 


(54) 


(t-i)pr(?)/’r(^),. 

- - A,  „ 

Di,„ 

(55) 

(x-i)pr(l)/'r(4)„. 

D,.  ®'- 

(56) 

a' 

-  I-,  ^l,m 

(57) 

B‘ 

(58) 
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AXIAL  FIELD  AND  GRADIENT  AT  THE 
CENTER  OF  THE  SPHEROID 


The  performance  of  the  permeable  spheroid  as  a  magnetic  sensor  core  is  determined  by  the 
combined  effects  of  the  permeability  and  the  demagnetizing  effect  of  the  body.  These  can  be 
examined,  for  both  the  field  and  gradient,  by  evaluating  the  two  expressions  at  the  origin.  Since 
the  expansion  for  <1>‘  is  identical  in  form  to  that  for  d>‘,  (Eq.  41)  can  be  written  in  terms  of 
cartesian  coordinates  with  the  aid  of  (Eq.  23-30)  used  in  reverse.  This  gives 


a;,,  Bl,  a;.o  9Ai, 

<P  =^— jc  +— -I — —2  + — —xz  + — —yz  + 


18512 


2  —  ■  2  ■  2 

a  a  a  a 


a  a 


(59) 


Then 


ad>‘ 

dXj 


(60) 


and 


a^d)‘ 

dXjdXi 


(61) 


Thus,  the  axial  field  is 


A;,o  9A1,  951,  48.41o 


-X  — 


■y  — 


(62) 


and  the  axial  gradient  is 


48Alo 


(63) 


Applying  (Eq.  35)  and  (Eq.  57),  the  ratio  of  the  interior  axial  flux  density  at  the  origin  to 
the  applied  flux  density  is  given  by 
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B‘(0.0,0)  H',.0  ^ 

- =  'l-r— = /X 

fis'  ^1.0 


and  the  ratio  of  the  axial  induction  gradients  everywhere  in  the  interior  is  given  by 


IIG33  ^2^0 

=  x^r— =  gx. 


II0G33  £^2.0 

The  Wronskian  Wi  „(^)  has  the  explicit  form 


( l+m+2'] 

<  2  J 

ri 

'  l  +  m  +  \  ' 

1  2  j 

(?-i)  rl 

.  2  /I 

1  2  j 

which,  for  m  =  0  reduces  to 


With  these  results,  the  quantities  /  and  g  are  given  by 

1 


/=- 


(x-1) 


[^(?-l)lnV|5T-?J  +  '' 


and 


(64) 


(65) 


(66) 


(67) 


(68) 


4-h3(t-l)^(3f^-l)lnVf{-6^]’ 

If  D  and  L  are  the  spheroid  diameter  and  length,  respectively,  then,  with  e  =  OIL, 

\  =  l/'Vl  - e^.  The  functions / and  g  are  plotted  in  Figure  1 .,  for  a  value  of  x  =  100.  It  is  clear 
from  these  curves  that  demagnetization  is  more  significant  for  the  gradient  than  for  the  field,  and 
an  LID  of  27  is  required  to  get  50%  of  the  benefit  of  the  higher  permeability. 
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